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INTRODUCTION 

Since  1970  eleven  stations  with  high-gein  (maguifice 
tions  up  to  150,000)  end  long  period  (peak  between  30  end 
50  seconds)  systems  recording  in  analog  and  digital  form 
have  been  put  into  operation  (Sevino  et  el.,  1972).  These 
stations,  that  will  be  referred  to  as  HGLP,  have  bean  used 
in  the  determination  of  source  mechanism  parameters,  compa 
ring  their  performance  with  that  of  the  World  Wide  Standard 
Station  Net  (WWSSN^  stations.  Locations  of  the  HGLP  stati- 
ons whose  data  have  been  used  in  this  study  ere  givon  in  te 
ble  1. 

A very  important  step  in  the  analysis  of  the  HGLP  da- 
ta is  the  preliminary  processing  which  involves  enalog  to 
digital  conversion,  filtering,  spectral  analysis  end  reduc- 
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tion  to  ground  motion.  The  HGLP  instruments  have  two  out- 
puts, one  enelog  (HGLP-ANL)  with  conventionel  photogrephic 
recording  the  other  digital  (HGPl-DIG)  on  magnetic  tape. 
Digitalization  of  the  analog  records  have  been  done  from 
enlarged  copies  of  the  microfilms.  In  thin  respect  the  da 
ta  from  earthquakes  relatively  large  end  near  to  the  sta- 
tion ere  practically  impossible  to  analyze.  Amplitudes  ere 
too  large,  lines  feinted,  and  traces  entengled.  The  digi- 
tal outputs  should  solve  this  problem,  but  it  creates  some 
of  its  own. 

In  order  to  study  the  capebilities  of  the  digital  out 
puts  of  the  HGLP  stations,  the  digital  tapes  in  file  at  the 
Lemont-Doherty  Observatory  were  studied.  Some  of  the  pro- 
blems found  were  the  following;  records  not  in  the  Lamont 
fils,  copies  defective,  failures  present  in  the  digital  sys 
tern,  and  too  loui  magni  ficetions.  The  process  of  finding 
the  proper  place  in  the  tapes  is  also  very  tima  consuming, 
since  there  la  not  a fast  search  system.  Another  problem 
encountered  in  using  HGLP-DIG  deta  is  the  varietion  of  the 
conversion  fector  to  ground  motion.  In  the  ALQ  (Albuquer- 
que Observatory)  station  the  values  found  in  10  calibrati- 
ons fror  Oct  73  to  Mar  74  give  for  the  magnification  et  30 
sac  in  counts  per  micron  values  ranging  from  802  to  1232  in 
the  Z component,  512  to  B52  in  NS  and  593  to  965  in  EtU  (A. 
Murphy , personal  comunication) . It  is  not  known  whether  the 
henges  were  the  result  of  modifications  made  to  the  system 
or  spontaneous  chenges.  If  the  letter  is  true,  large  errors 
can  be  introduced  when  using  digital  deta  from  stations  whe 
re  calibration  is  not  repeated  so  often.  This  lack  of  sta- 
bility in  the  digital  system  would  be  reflected  in  the  compa 
rison  of  the  source  parameters  which  require  ground  motion 
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reduction  as  determined  from  analog  and  digital  records. 

SOURCE  PARAMETERS  AN 3 THE  SPECTRUIT)  OF  SEISMIC  WAVES 

For  many  yeers  all  that  could  be  found  about  the  me- 
chanism of  earthquakes  was  the  orientation  of  the  two  nodal 
planes  of  P waves,  one  of  them  corresponding  to  the  plane  of 
fracture  and  the  energy  released  as  seismic  waves  derived 
usually  from  magnitude.  Progress  in  the  study  of  focal  me- 
chanism in  the  last  ten  years  has  led  to  the  definition  and 
determination  of  several  parameters  describing  more  fully 
the  processes  at  the  source  of  earthquakes,  euch  as  the 
fracture  dimensione,  dislocation,  seismic  moment,  average 
stress  and  stress  drop.  Dislocation  dynamic  models  of  the 
source  of  earthquakes  predict  the  shape  of  the  spectrum  of 
seismic  waves  for  the  far  displacement  field.  From  the  pa- 
rameters which  define  the  form  of  such  spectrum  certain  im- 
oortant  characteristics  of  the  dimensions  and  stress  condi- 
tions at  the  source  can  be  determined.  many  authors  have 
in  the  recent  years  treated  this  problem,  among  them,  Aki 
(1967,  1972)  Ben  menahen  (1961,  1965),  Brune  (1970),  Hanks 
and  Thatcher  (1972),  Hanks  and  Ufyss  (1972),  Dahlen  (1974). 

In  general  lines  the  spectrum  of  seismic  waves  at  the 
^ar  pield  can  be  described  as  consisting  of  a low  frequency 
flat  part  of  amplitude  00  followed  by  a decay  of  amplitude 
with  frequency  starting  at  a particular  frequency  f called 
the  "corner  frequency".  The  way  epectral  amplitude  decay 
with  frequency  is  governed  by  the  parameter  £ . For  a par 
ticular  case  when  6=1,  the  spectral  amplitudes  decay 
approximately  as  the  inverse  of  the  square  of  frecuency. 

This  case  represents  a dislocation  in  which  the  strees  drop 
is  complete.  These  three  parameters  of  the  spectrum  can  be 


4 


TABLE  1 


High-Gain  Long-Period  station  locations. 


Station 

Code 

Station  Name 

Location 

OGD 

Ogdensburg 
New  Jersey 

4l°04'0ON  ; 74°J7'00W 

CHG 

Chiang  Mai 
Thailand 

1B°47'24"N;  98°58'37"W 

EIL 

Eilat 

Israel 

29°33,00"N;  34°57,00”E 

FBK 

Fairbanks 

Alaska 

64°53 ' 58 ' 'N;  l4a°00'20"W 

TLO 

Toledo 

Spain 

39°51'36"N;  4°00'53,,l* 

KON 

Kongsberg 

Norway 

59°38'57"N;  9°37’55"E 

ALQ 

Albuquerque 
New  Mexico 

34°56 ' 30  " N ; 106°27 ' 30 ’ ' W 

Elevation 

Meters 

-373 

416 

330 

465 

200 

1053 
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related  to  source  characteristics.  The  spectral  amplitude 
level  al  low  frequencies  Qo  is  related  to  the  seismic  mo- 
ment, uihich  for  i fracture  is  defined  as  the  product  of  the 
rigidity,  area  and  average  dislocation.  The  corner  frequen 
cy  fQ  is  related  to  the  dimensions  of  the  feult,  and  e.  with 
the  retio  of  the  stress  drop  to  the  effective  stress.  In 
principle  spectre  from  either  budy  or  surface  weves  can  be 
used  in  the  determination  of  source  parameters.  Howaver, 
surface  weves  are  preferred  in  the  determination  of  seismic 
moments  of  large  shallow  shocks  since  for  these  cases  the 
long  period  soectrel  date  of  body  waves  are  mixed  with  la- 
ter errivels.  For  smaller  and  deeper  events  body  waves  are 
found  more  useful.  For  the  determination  of  source  dimen- 
sions from  the  corner  frequency  is  preferable  to  use  body 
waves  since  fQ  is  generally  et  a frequency  where  surface  wa 
ves  depend  greatly  on  propagation  conditions  (Hanks  and  U/yss 
1972).  For  the  magnitude  range  (5  <m  ^ 6)  of  the  shocks 
of  this  stud/,  and  shallow  depths,  surface  and  body  waves 
can  be  equally  used  for  the  determination  of  seismic  mo- 
ments. Source  dimensions  cen  be  calculated  from  the  corner 
frequency  of  shear  waves  ana  also  from  surface  waves,  but 
not  fro.r  the  corner  frequency,  but  from  the  directivity  func 
tion.  In  this  way  for  the  same  shock  source  parameters  can 
be  determined  independently  from  body  and  surface  waves. 

Source  parameters  from  spectra  of  Rayleigh  waves 

Following  the  work  of  Ben-fflenehen,  Rosenman  and  Har- 
krider  (1970)  the  seismic  moment  M0(R)  cen  be  expressed  in 
terms  of  the  spectial  amplitude  et  the  low  frequencies  of 
the  vertical  component  of  the  Reyleigh  wave  £2*00  in  the 
f orm; 
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M0(fO  = 


_ HjlM) 


(1) 


V sin  .A_ 
a. 


(^r*r  + srjr"  1 j 


u)hBre  K is  the  rigidity  modulus,  the  attenuation  coe 

fficient  of  Rayleigh  waves,  A the  distance,  a the  earth 

radius,  and  PR,  pR,  SR,  sr,  Qr,  qR,  are  the  transfer  func- 

tions for  a stratafied  earth  corresoonding  to  a double  cou- 
ple source.  The  values  of  the  parameters  of  this  equation 
are  given  in  convenient  table  form  by  Ben-menahen,  Rosen.nan 
and  Harkrider  (1970), 

The  fault  length  can  be  obtained  assuming  a propaga- 
ting fracture  of  rectangular  form,  from  the  "directivity 
function".  This  function  has  been  defined  by  Ben-menahen 
,1961)  as  the  ratio  of  the  spectral  amplitudes  of  surface 

waves  corresponding  to  rays  leaving  the  source  in  opposite 

directions.  The  limitation  of  using  only  data  from  stations 
190  degrees  apart  in  azimuth  is  not  necessary  for  faults 
with  large  strike  slip  component  (Udfas,  1971).  If  a is 
the  angle  between  the  azimuths  of  two  stations,  the  length 
of  the  fault  b is  given  by 


t = 


~ ~ ("os  (fy  ♦-  ) 


°,here  fmax  and  f min  are  the  fre9uencies  of  the  first  maxi- 
mum or  minimum  of  the  directivity  function,  v the  rupture 

velocity,  cr  the  ! yleigh  wave  phase  velocity  and  9 the 
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azimuth  from  the  fault  plane  to  one  of  the  stations. 

In  order  to  compare  the  results  of  source  parameters 
determined  from  the  spectrum  of  Rayleigh  waves  with  magni- 
tude, is  convenient  to  calculate  it  also  from  surface  waves. 
The  following  formula  has  been  used  hare 

Ms  = -1.66  too  A + 3 •‘3  (3) 

where  A is  the  amplitude  in  microns  of  surface  waves  corres- 
ponding to  20  s period,  T the  period  and  the  epicentral 
distance  in  degrees.  From  this  value  of  magnitude,  the  ener 
gy  can  be  calculated  using  the  relation  (Richter,  1958) 

leg  E5  = H.&  + 1.5  t^s 

Source  parameters  from  the  spectra  of  5 waves 

Accepting  the  general  validity  of  Brune's  model  (Bru- 

ne,  1970  and  1971)  and  the  physical  interpretation  given  by 

him  to  the  three  parameters  and  £ , derived  f rom 

the  spectra  of  shear  waves,  they  can  be  used  to  determine 

the  seismic  moment  m and  the  dimension  r of  the  equivalent 

o 

circular  dislocation.  Here  £ is  taken  to  be  unity,  assu- 
ming that  the  stress  drop  is  complete. 

The  seismic  moment  m (s)  can  be  determined  from  the 

o 

spectral  amplitude  level  at  low  frequencies  Qt>(s')  of  the 

shear  wives  using  the  relation  (Hanks  and  UJyss,  1972) 

IS  R 

Mo(S)  = * k yp3  R £2o(s)  k £ (5) 

Where  fb  is  the  shear  wave  velocity  near  the  sourca,  and  R 
the  hypocentral  distance.  The  attenuation  coefficient  •g  is 
given  by 

ii 

* = Tfyp 
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o/here  T is  the  period  and  Qp,  the  average  quality  factor 
for  shear  waves.  The  factor  k corrects  for  the  effect  of 
the  radiation  pattern  and  the  reflection  at  the  free  surfa 
ce  and  is  given  by 

k = - 

where  Is  the  normalized  radiation  pattern  for  the 

S - waves  from  a double  couple  source. 

The  corner  frequency  of  the  shear  waves  spectrum 
fQ(S)  is  related  to  the  radius  of  the  equivalent  circular 
dislocation  r by 


w 

For  a comparison  between  the  dimensions  defined  by  r and  b, 
one  may  consider  that  2r  represents  thB  width  of  the  rectan 
gular  fault  and  '.his  as  1/3  of  the  length  b.  Another  rela- 
tion can  be  derived  from  equal  area  considerations 


which  results  in 


V ~ k = O.  525  t 


The  energy  radiated  as  seismic  waves  Es  is  given  as  a 
fraction  of  the  total  energy  E by 

Es  * 1 E (9) 

where  ^ is  the  seismic  efficiency  factor.  This  enerqy  f 

s 

can  be  calculated  from  the  body  wave  magnitude  by  the 
standard  relation  (Richter,  1958) 


log  Es  - 5.8  + 2.5 


(10) 
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In  terms  of  flQ(S)  and  fQ(S)  the  energy  is  also  given  appro 
ximately  by  (UJyss  and  Hanks,  1972) 


Es  SL  8.53  KZf  (b  R2Qo(S)|0(s'v  k2 


(ID 


Average  dislocation,  average  stresi  and  stress  drop 


From  the  analysis  of  either  surface  or  body  uiavos  we 
can  independently  determine  the  seismic  moment,  fault  dimen- 
sions (radius  or  length)  and  seismic  energy.  From  these  pa- 
rameters others  nan  be  derived  that  are  useful  in  order  to 
describe  the  nature  of  a seismic  source. 


For  the  source  of  an  earthquake  considered  as  a sudden 
n location  on  a plane  surface,  the  seismic  moment  is  defi- 
ned as  the  product  of  the  area  of  the  fault  A by  the  rigidi 
ty  )*  and  the  average  dislocation  or  displacement  over  the 
fault  surface  w (Wyss  and  Brune,  I960) 


= ,u  /\  w 


(12) 


From  this  relation  the  average  dislocation  is  given  by 


\KJ  - 


Mo 


A A 


(13) 


and  for  the  rectangular  fault  model  and  the  circular  one 
used  in  the  surface  and  body  wave  analysis  one  can  write 


w(fO  - 5 Mo(R) 


(14) 


w (5)  - 

^ tc  r2- 


(15) 


The  difference  in  shear  stress  acting  on  a fault  surface  be 
fore  and  after  the  occurrence  of  an  earthquake  is  generally 


lo 


referred  to  as  the  "stress  drop"  A <J=  cr<-crz  This  is  rela- 
ted to  the  maximum  dislocation  w^  through  the  equation  (Bru 
ne  and  Allen,  1967) 

Act  - ^ Ms  if-  (16) 

d 


u/here  d is  the  width  of  the  fault  and  ^ a constant  ranging 
from  1/2  to  2tc/3,  depending  on  fault  geometry.  For  earth- 
quakes with  large  strike-clip  component  of  motion  ^ can  be 
taken  as  1/2.  The  maximum  dislocation  on  a fault  is  usually 
taken  as  4/3  the  average  value  w . Taking  these  relations 
and  the  expression  for  ui  from  equation  (13)  ujb  obtain  for 
the  stress  drop  in  a vertical  rectangular  3trike-slip  fault 
with  width  equal  to  1/3  its  length 


ACT  (R)  = 


' k5 


(17) 


This  relation  is  similar  to  the  one  given  by  Brune  (1971) 
for  a circular  dislocation 

Acr(s^i  = 7 H.C5) 

16  rs  (IB) 

As  the  stress  drop  AH  ) depends  on  the  square  of  the 
fault  width,  it  is  greatly  influenced  by  the  assumed  rela- 
tion between  fault  length  and  width. 

The  total  energy  E released  by  slipping  on  a fault 
can  be  written  as  the  product  of  the  average  shear  stress 
a = ( CT,  -h  CTZ ) acting  on  the  fault  plane  by  the  area  and 

the  average  dislocation 

E = CT  A \AJ 


(19) 
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I 

From  this  relation  the  "apparent  average  stress"  rjff  , cr 
product  of  the  average  stress  times  the  seismic  efficiency 
factor  can  be  obtained  as, 

rp  - (20) 

Since  Eb  end  fnQ  can  be  independently  calculated  from  body 
and  surface  waves,  me  can  also  calculate  two  values  of  the 
apparent  stress.  The  interpretation  of  the  ratio  of  the  seis 
mic  energy  to  the  moment  depends  on  the  source  model.  Ran- 
dall (1972)  interprets  it  as  being  related  to  the  stress  drop 
instead  of  the  average  stress. 

SOURCE  PARAMETERS  OF  FOUR  EARTHQUAKES 

Four  earthquakes  have  been  selected  in  the  same  magni- 
tude range,  5 < Itl  < 6 ; dates,  epicenters  and  origin  times 
are  given  in  table  2 . Three  are  located  in  the  mid-Atlan- 
tic ridge  and  one  in  the  Aleutian  islands.  Focal  mechanism 
using  first  motions  of  P wave  have  been  determined;  solu- 
tions are  shown  in  figures  1,  2,  3,  and  4,  and  orientation 
(strike  and  dip)  of  the  two  planes  given  in  table  2.  Solu- 
tion for  shock  AZOR  represents  normal  faulting  along  planes 
oriented  approximately  north-south.  This  solution  is  consist 
tent  with  the  type  of  faults  present  in  an  expanding  ridge. 

The  solution  for  ALEU  is  not  well  defined.  The  data  allows 
two  possible  solutions  one  involving  strike-slip  motion  on 
vertical  planes  striking  north-south  and  east-west  and  the 
other  reverse  faulting.  For  the  first  type  of  mechanism  the 
north-south  plane  which  has  been  taken  as  the  fault  plane  is 
consistent  with  an  "arc-arc"  transform  fault.  This  type  of 
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faulting  is  common  in  the  Aleutian  Islands  (stauder,  1968). 
The  second  solution  has  no  easy  tectonic  interpretation.  Re 
verse  faulting  in  the  Aleutian  region  give  planes  with  an 
east-west  orientation  and  underthrusting  of  the  Pacific  block 
under  the  islands  consistent  with  the  subduction  movement 
of  the  oceanic  plate.  The  solution  for  NAT  1 and  NAT2  repre 
sent  right-lateral  strike-slip  movement  on  an  east-west 
fault,  motion  wich  is  consistent  with  transform  faults  of 
the  mid-Atlantic  ridge.  In  conclusion  except  for  AZOR  the 
mechanism  have  large  strike-slip  component  of  motion. 

In  figure  1,  2,  3 and  4 the  position  of  the  stations 
used  in  the  Rayleigh  wave  analysis  are  shown  on  the  rim  of 
the  projection  and  that  of  the  stations  used  in  the  body 
wave  analysis  with  a large  open  circle. 

Determination  of  seismic  moment  and  dimensions 
Spectra  of  Rayleigh  waves 

Spectra  of  Rayleigh  waves  have  been  calculated  for  a 
tona^  of  21  stations,  of  these  6 HGLP  and  15  IMUSSN.  Exam- 
ples of  the  spectra  are  given  in  figures  5,  6,  and  7.  All 
spectra  show  a well  defined  constant  level  of  spectral  am- 
plitude for  low  frequencies.  Figure  9 shows  the  soec 

tra  of  the  Rayleigh  waves  in  ALQ  from  the  records  of  the 
UlUJSSN  and  HGPL  (analog  and  digital).  The  two  analog  records 
give  very  similar  spectra. 

Determinations  of  the  seismic  moment  ff]Q  has  been  done 
using  equation  (l)  and  results  are  given  in  table  3.  Va- 
lues are  relatively  stable  and  show  a good  agreement  bet- 
jeen  the  determinations  based  on  data  from  tl/UJSSN  stations 
and  HGLP  stations.  The  spectral  amplitudes  have  been  mea 
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TABLE  3 


Spectral  amplitudes  of  Rayleigh  waves  and  seismic  moments. 

T - 100  S 


T - 50  S 


STATION 

(km) 

M^IO^4  dyne-cm) 

M (I0f>4  d 

ov 

AZOR 

AAE-WWSSN 

8025 

O.G 

1.4 

0.6 

2.9 

8UL 

9325 

1.8 

6.4 

2.4 

38.0 

COP 

4290 

2.4 

3.4 

2.6 

3.0 

LON 

7006 

1.3 

2.1 

1.5 

5.3 

NUR 

5079 

3.1 

4.4 

4.4 

6.1 

ALEU 

ALQ-KWSSN 

6057 

3.0 

2.9 

1.4 

0.6 

MAT 

3602 

0.8 

0.8 

0.3 

0.1 

NOR 

5213 

0.2 

0.2 

0.03 

0.01 

NAT1 

EIL-HGLP(AN) 

5706 

C.4 

0.6 

0.16 

0.2 

KON 

2293 

1.0 

1.6 

0.19 

0.2 

OGD 

3758 

2.3 

3.1 

0.45 

0.3 

GDH-WWSSN 

3054 

6.0 

8.3 

2.3 

1.3 

KON 

2993 

2.0 

3.3 

1.8 

1.5 

VAL 

1530 

5.5 

8.0 

3.3 

1.4 

NAT  2 

ALQ-HGLP(AN  ) 

6041 

2.9 

9.0 

1.4 

1.6 

-ttGLP  ( DTS ) 

1.1 

3.4 

0.2 

0.2 

-WWSSN 

9.2 

29.0 

3.2 

3.6 

kon-hglp(an  ) 

4671 

1.3 

1.9 

0.06 

0.03 

TLO 

3285 

2.9 

4.4 

0.6 

0.6 

LP8-WWSSN 

6263 

0.3 

0.8 

0.05 

0.2 

15 
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sured  at  two  periods,  100  and  50  seconds.  (Tlomerts  based 
on  the  amplitudes  at  50  s period  are  lower  on  the  average 
than  those  of  the  100  s.  Taken  average  values  the  lowest 
moment  is  that  of  NA1 1,  but  the  difference  is  so  small 
that  practically  all  moments  car,  be  say  to  be  equal. 

To  examine  the  scatter  of  the  data  and  the  influence 
of  the  radiation  pattern  factor  we  have  plotted  in  figure 
9 the  values  of  the  spectral  amplitudes  corrected  by  sprea 
ding  and  attenuation  versus  distance.  In  the  figure  the 
curves  corresponding  to  constant  moment  are  shown.  These 
have  been  calculated  for  a theoretical  model  with  vertical 
strike-slip  faulting  at  a surface  focus  (30  km)  and  corres 
ponds  to  soectral  amplitudes  recorded  a*-  45°  from  the  fault 
plane.  Scatter  of  the  observed  spectral  amplitudes  amounts 
on  the  average  to  one  order  of  magnitude.  Spectral  ampli- 
tudes at  50  s are  more  scattered  than  those  at  100  s. 


If  moments  are  calculated  from  the  spectra]  amplitu- 
des at  100  s using  the  diagram  of  figure  9 prescinding  of 
the  focal  mechanism,  the  average  for  each  shock  are  in  good 
agreement  with  the  average  values  obtained  from  table  3 who 


factors 

have  been  considered 

in  th->  calculations 

shock 

^(calculated) 

(^(diagram) 

24 

10  dyne-cm 

24 

10  dyne-cm 

AZ0R 

3.5 

2.4 

ALEU 

1.3 

2.4 

NAT  1 

4.1 

1.9 

NAT  2 

3.9 

3.3 

Differences  are  within  the  range  of  error  usually  admitted 
for  this  parameter.  Provided  that  stations  are  not  near  a 
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nodal  plane  and  that  mechanisms  hove  large  strike-slip  com- 
ponont  of  motion  curves  of  figure  9 should  give  on  the  ave- 
rage sufficiently  good  estimates  of  the  seismic  moment. 

Fault  lengths  have  been  calculated  for  the  four  shocks 
from  tho  directivity  function  according  to  equation  (2).  The 
rupture  velocity  has  been  taken  as  1.5  k:n/sec.  An  example 
of  the  directivity  function  for  NAT2  calculated  from  the 
spectra  of  HGLP  stations  of  KON  and  TLO  is  shown  in  figure 
10.  Values  of  b are  given  in  table  4.  The  rupture  veloci 
ty  has  been  chosen  to  keep  the  fault  dimensions  in  agreement 
with  the  results  from  body  waves.  However,  increasing  v to 
2 km/sec,  a more  commonly  accepted  value,  increases  b in  on- 
ly about  10  km. 

Spectra  of  Shear  UJaves 

Analysis  of  body  waves  has  been  limited  to  only  two 
shocks  and  the  study  of  S waves.  Only  for  ALEU  and  NAT2  the 
S wave  pulses  were  sufficiently  clear.  Spectra  are  shown  in 
figures  11,  12,  13  and  14.  Data  from  HGLP  stations  have  Deen 
ffund  very  useful,  since  long  period  5 waves  are  there  usua- 
lly well  recorded.  The  spectral  level  at  low  frequencies 
are  in  general  well  defined.  Corner  frequencies  are  loss  re 
liable,  with  a certain  amount  of  uncertainty  with  resoect  tn 
the  position  of  the  point  of  intersection  of  the  constant  le 
vel  aith  the  sloping  line  of  amplitude  with  frequency. 

Values  for  seismic  moment  and  radius  of  the  circular 
fracture  have  been  calculated  using  equations  (5)  and  (6). 
Results  are  shown  in  table  5.  moments  are  similar  to  those 
found  from  Rayleigh  waves  spectra,  though  somewhat  lower.  Fi 
gure  15  shows  the  spectral  amplitudes  distance  relation  with 


TABLE  4 

Length  of  fault  from  diractivity  of  Rayleigh  waves. 


Shock 

Stations 

b (Km) 

AZOfi 

BUL  / LON  (WWSSN) 

27 

ALEU 

ALQ  / NOR  (WWSoN) 

54 

N 

47 

NAT1 

AQU  / BEC  (WWSSN) 

43 

NAT  2 

KON  / TLO  (HGLP  ) 

31 

KON  / ALQ  (HGLP  ) 

(fR) 

TLO  / ALQ  (HGLP  ) 

42 

TABLE  5 


Seismic  moments  and  diameter  of  fault  from  shear  waves. 


STATION 

fl  (km) 

of-®) 

24 

M^lO  dyne-cmj 

f0(cPs) 

2r(km) 

ALEU 

ALQ-HGLP(  ANAL  )— ■(  03  ) 

6050 

0.31 

0.2 

0.092 

15.4 

F0K-HGLP  ( AN  AL  )-(03) 

2415 

- 

- 

0.100 

14.2 

-(02) 

•1 

- 

- 

0.046 

31.0 

CGH-HGLP(DTS  )-(A3) 

7050 

0.3 

0.0 

0.039 

36.3 

COL-WWSSN  “(05) 

2427 

2.7 

0.3 

0.051 

27.7 

-(06) 

*1 

1.5 

0.002 

17.3 

COR-WWSSN  -(06) 

4266 

1.3 

0.2 

0.051 

27.7 

DUG-WWSSN  “(06) 

5253 

0.26 

0.2 

0.002 

17.3 

NAT  2 

ALQ-HGLP(DTS  )-(A2) 

6041 

1.62 

1.6 

0.070 

13.0 

K0N-HGLP(ANAL)-(05) 

4672 

0.57 

0.4 

0.002 

17.3 
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curves  of  constant  moment.  Curves  have  been  calculated  with 
out  the  radiation  pattern  term,  taking  instead  the  mean  va- 
lue of  the  radiation.  The  spread  of  the  observed  spectral 
amplitudes  is  much  less  than  for  Rayleigh  wave  data  (figure 
9).  The  correction  for  the  radiation  pattern  is  in  this  ca- 
se not  very  important  and  values  of  moments  can  be  calcula- 
ted with  sufficient  approximation  without  that  term.  Compa- 
ring the  diameter  obtained  for  the  circular  crack  with  the 
length  of  the  rectangular  fault  used  in  the  directivity  cal- 
culations, it  corresponds  to  the  equivalent  dimension  for 
equal  area  correspondance  taking  a fault  width  of  2/3  the 
fault  length. 


Seismic  energy,  stress  drop  and  average  stress 


To  complete  the  determination  of  source  parameters  the 
radiated  seismic  energy  has  been  calculated  using  the  equa- 
tions (A)  and  (10)  from  the  values  of  magnitude.  Surface  w£ 
ve  magnitudes  have  been  calculated  from  amplitude  readings 
at  several  stations  and  by  means  of  equation  (3).  Results 
are  shown  in  table  6.  For  body  wave  magnitude  that  given  ty 
USCGS  has  been  adooted.  Energies  calculated  from  m^  are 

lower  than  those  from  ffl  . This  is  due  to  the  fact  that  mL 

s b 

is  lower  than  the  body  wave  magnitude  corresponding  to 
our  values  of  Iti  according  to  the  relation 

m()  - 2.5  * 0.63  (21) 


Stress  drops  AO'  have  been  calculated  from  Rayleigh  w?a 
ve  data  using  equation  (17)  and  from  S wave  data  using  (18). 
Results  are  shown  in  table  6.  Stress  drops  are  rather  low, 
specially  for  ALEU  (R)  due  to  the  high  value  obtained  for  b. 
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Source  parameters. 


TABLE  6 


E 

s 

M 

o 

b,2r 

A 

AC 

w 

SHOCK 

M 

s’  b 

.19 

, _24 

2 

10  ergs 

10  dyne— cm 

km 

km 

bars 

bars 

cm 

Rayleigh 

waves 

AZOR 

5.5 

11 

3.5 

27 

243 

1,1 

10.4 

4.3 

ALEU 

5.5 

11 

1.3 

51 

B67 

.06 

27.9 

0.4 

NATl 

5.0 

2 

4.1 

43 

616 

0.3 

1.6 

2.0 

NAT  2 

5.4 

9 

3.9 

36 

432 

0.5 

7.6 

2.7 

S - waves 

ALEU 

5.5 

1 

0.3 

23 

416 

0.1 

11 

0.2 

NAT2 

5.5 

1 

1 

15 

176 

1.2 

3.3 

2.0 

- ' 
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In  figure  16  seismic  moments  have  been  plotted  versus  fault 
area.  Curves  of  constant  stress  drop  have  been  also  drawn 
for  Brune's  circular  crack  model.  The  stress  drops  found 
here  vary  between  0.1  and  1 bar  and  are  lower  by  an  order 
of  magnitude  than  tne  expected  from  Brune's  model.  In  the 
same  figure  me  have  plotted  values  determined  by  other  au- 
thors using  the  same  methods  as  here  (300  table  7).  uje 
find  in  these  values  again  8 good  agreement  between  calcu- 
lations using  Rayleigh  and  S waves.  Stress  drops  are  bet- 
ween 0.1  and  10  bars  and  also  a little  lower  than  the  ex- 
pected from  the  model. 

The  ratios  of  the  seismic  energy  to  mement  have  been 
interpreted  as  giving  a measurement  of  the  apparent  average 
shear  stress  on  the  fault  zone.  The  values  obtained  here 
vary  from  1 to  30  bars.  Since  ^ must  be  less  than  unity, 
these  are  minimum  estimates  of  the  average  stress.  The  or- 
der of  magnitude  obtained  agrees  with  that  generally  accep- 
ted for  the  strength  of  the  earth  crust  under  horizontal 
shear  stress  (Chinnery,  1964).  Randall's  contention  that 
this  ratio  actually  represents  1/2  ACT  ^ does  not  t,gree 
with  our  results  where  its  values  are  higher  than  those  ob- 
tained for  A'J  from  the  moment  and  dimensions  of  the  fault. 
Low  stresses  are  expected  at  the  zone  of  oceanic  ridges, 
which  will  explain  the  values  found  for  AZOR,  NAT1  and  NAT2. 
The  stress  for  ALEU  is  higher,  but  due  to  the  uncertainties 
intrinsic  to  this  type  of  determinations  it  is  hard  to  say 
if  this  difference  is  significant. 

Accepting  the  interpretation  given  to  the  values  of 
the  stress  drop  and  average  stress,  the  relation  between 
them  will  give  a fractional  s^ess  drop  of  10^  for  seismic 


Source  parameters  from  other  authors 
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efficiency  equal  to  unity.  A similar  result  was  found  by 
UJyss  and  Brune  (1958)  for  earthquakes  in  the  California  Ne- 
vada region. 

Average  dislocations  have  been  determined  from  equa- 
tion (14)  and  (15).  Values  range  from  0.2  to  4 cm.  These 
values  agree  with  those  found  by  other  authors  for  shocks 
of  similar  magnitude. 

CONCLUSirNS 

Records  from  high-gain  long-period  (HGLP)  stations 
have  been  found  useful  in  the  analysis  of  source  parameters 
of  earthquakes  of  magnitude  about  5.  For  this  magnitude 


level  anJ  distances  between  2000  and  8000  km  Rayleigh  wa- 
ves are  well  recorded  and  in  many  cases  clear  S wave  pulses 
can  be  found.  These  data  is  well  suited  for  determination 
of  source  parameters  from  spectra  of  both  types  of  waves. 


Independent  determinations  of  source  parameters  from 
Rayleigh  and  shear  wave  data  show  good  agreement.  In  each 
case  one  can  find  independent  estimates  of  seismic  energy, 
moment  and  dimensions  of  the  fault  and  from  these  derive 
values  for  stress  drop,  apparent  average  stress  and  average 
dislocation . 


For  surface  shocks  of  magnitude  between  5 and  6 u»e 
have  found  consistently  low  stress  drops  between  0.1  and 
1 bars  which  compared  with  the  apparent  average  stress  gi- 
ve maximu*  fractional  drops  of  stress  of  10$.  Faults 

2 

areas  are  between  200  and  800  km  and  total  seismic  moments 
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of  the  order  of  10  dynes-cm. 
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1.  Focal  mechanism  solution  for  tha  Azores  earthquake  of  18 
Now  1970  (AZOR).  Shown  the  positions  of  stations  used 
in  Rayleigh  wave  analysis. 

2.  Focal  mechanism  solution  for  the  Aleutian  Islands  earth- 
quake of  3 Jan  1972  (ALEU).  Shown  the  positions  of  the 
stations  used  in  the  Rayleigh  wave  (rim)  and  shaar  wave 
(inside)  analysis. 

3.  Focal  mechanism  solution  for  the  North  Atlantic  earthqua 
ke  of  6 Jun  1972  (NAT1).  Shown  the  position  of  the  eta” 
tions  used  in  Rayleigh  wave  analysis. 

4.  Focal  mechanism  solution  for  the  North  Atlantic  earthqua 
ke  of  6 Jun  1972  (NAT2).  Shown  the  position  of  the  sta- 
tions used  in  the  Rayleigh  wave  (rim)  and  shear  wave  (in 
side)  analysis. 

5.  Spectrum  of  Rayleigh  wave  at  UMUSSN  station  of  BUL  for  the 
18  Nov  1970  (AZOR)  earthquake. 

6.  Spectrum  of  Rayleigh  wave  at  the  UJtUSSN  station  of  NOR 
for  the  3 Jan  1972  (ALEU)  earthquake. 

7.  Spectrum  of  Rayleigh  wave  at  the  HGLP  station  of  KON  for 
the  13  Way  1972  (NAT1)  earthquake. 

8.  Spectra  of  Rayleigh  wave  for  the  HGLP  (digital  and  analog 
records)  and  UJtUSSN  stations  at  ALQ  for  the  6 Jun  1972 
(NAT2)  earthquake. 

9.  Relation  between  the  spectral  amplitudes  corrected  for  at 
tenuation  and  spreading  and  distance  for  Rayleigh  waves. 
Curves  of  constant  seismic  moment  given  in  dyne-cm. 

10.  Directivity  function  of  Rayleigh  waves  for  the  HGLP  sta- 
tions of  KON  and  TLO. 

11.  Spectrum  of  shear  waves  at  the  HGLP  station  of  F8K  for 
the  3 Jan  1972  (ALEU)  earthquake. 

12.  Spectrum  of  shear  waves  at  HGLP  station  of  ALQ  for  the  3 
Jan  1972  (ALEU)  earthquake. 


Spectrum  of  shear  waves  at  the  U1UISSN  station  of  COR  for 
the  3 Jan  1972  (ALEU)  earthquake. 

Spectrum  of  shear  waves  at  the  HGLP  (digital  data)  sta- 
tion of  ALQ,  for  the  6 June  1972  (NAT2)  earthquake. 

Relation  between  spectral  amplitudes  corrected  for  atte 
nuation  and  spreading  distance  for  shear  waves.  Curves 
of  constant  seismic  moment  are  given  in  dyne-cm. 

Relation  between  seismic  moment  and  fault  area.  Circles 
(open  Rayleigh  and  black  shear  wave  data)  from  table  6 
(this  study).  Triangles  (open  Rayleigh,  and  black  shear 
wave  data)  from  table  7.  Curves  of  constant  stress  drop 
in  bars. 
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18  NOV  1970 


Figure  !•  Focal  mechaniam  solution  for  the  Azores  earthquake  of  18 
Nov  1970  (AZOR).  Shown  the  positions  of  stations  used 
in  Rayleigh  wave  analysis. 


3 JAN  1972 


2.  Focal  mechanism  solution  for  the  Aleutian  Islands  earth 
quake  of  3 Jan  1972  (ALEU).  Shown  the  positions  of  the 
stations  used  in  the  Rayleigh  wave  (rim)  and  shear  wave 
(inside)  analysis. 
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Figure  3.  Focal  mechanism  solution  for  the  North  Atlantic  earthqua 
ke  of  6 Dun  1972  (NATl).  Shown  the  position  of  the  sta- 
tions used  in  Rayleigh  vave  anelysia. 


Figure  4.  Focal  mechanism  solution  for  the  North  Atlantic  earthqua 
ke  of  6 3un  1972  (NAT2).  Shown  the  position  of  the  sta- 
tions used  in  the  Rayleigh  wave  (rim)  and  shear  wave  (in 
side)  analysis.  “ 
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Figure  6.  Spectrum  of  Rayleigh  wave  at  the  IIAUSSN  station  of  NOR  for 
the  3 Jan  1972  (ALEU)  earthquake. 
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Figure  8.  Spectra  of  Rayleigh  wave  for  the  HGLP  (diqital  and  analog 

records)  and  U/UiSSN  stations  at  ALQ  for  the  6 3un  1972 
(NAT2)  earthquake. 
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Figure  14.  Spectrum  of  shear  eaves  at  the  HGLP  (digital  data)  sta- 
tion of  ALQ,  for  the  6 June  1972  (NAT2)  earthquake. 
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Figure  15.  Relation  between  spectral  amplitudes  corrected  for  atte 
nuation  and  spreading  distance  for  sheer  weves.  Curvps 
of  constant  seismic  moment  ere  given  in  dyne-cm. 
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Figure  16.  Relation  between  seismic  moment  and  fault  area.  Circles 
(open  Rayleigh  and  black  shear  wave  data)  from  table  6 
(this  study).  Triangles  (open  Rayleigh,  and  Black  shear 
wave  data)  from  table  7.  Curves  of  constant  stress  drop 
in  bars. 


1. 

ATTENUATION  OF  SEISMIC  WAVES  DETERMINED  FROM  HIGH-GAIN 
LONG-PERIOD  STATIONS 

Alfonso  L6pez-Arroyo  and  Julio  Mezcua 
Instituto  GeogrSf ico  y Catastral 

SUMMARY 


A comparison  is  made  of  attenuation  values  obtained  from  data 
recorded  by  High  Gain  Long  Period  (HGLP)  instruments  and  by 
stations  of  the  Standard  network.  Attenuation  refers  to  surfa- 
ce Rayleight  waves  generated  by  nuclear  explosions  and  the  fre- 
quency range  considered  extends  from  0.1  to  0.01  Hz.  For  the 
distances  involved  in  this  study  the  use  of  such  energy  sources 
eliminates  the  uncertainites  derived  from  lack  sufficent  know- 
ledge of  the  radiation  pattern,  as  this  can  be  assumed  to  be 
spnerlcally  symmetrical.  Attenuation  coefficient  for  every  fre- 
quency is  given  as  the  slope  of  a least  square  linear  fitting 
of  the  spectral  amplitud  vs  distance.  Values  obtained  correspond 
to  paths  from  Novaya  Zemblya  to  Eastern  USA,  Canada  and  Europe. 
Results  from  HGLP  show  large  scatter  which  reflect  the  insuffi- 
cient control  of  the  instrument  response  for  the  system.  Had 
this  problem  been  adequately  solved,  the  HGLP  data  would  have 
been  far  superior  than  standard  data  for  the  attenuation  study, 
due  to  their  larger  dynamic  range  and  readness  to  use  on  the 
computer. 


i 


i 


In  the  present  study,  we  are  concerned  with  the  decreases 
in  the  amplitude  density  spectrum  with  distance.  That  is, 
we  wish  to  know  how  the  amplitude  density  at  a given  period 
(or  frequency)  decreases  with  increasing  distance. 

For  the  fundamental  mode,  a seismic  surface  disturbances, 
can  be  written  as 

( _££iidL_  expf-r^l  e*p  f i («ot+0)]  dio  /n 
2 TT  J^v/S  ih  0 * ' 

where 

V/(x,w)  = «p[-K*] 

VSi*©- 

and  ^ is  a phase  lag.  The  quantity  V is  the  Fourier  trans- 
form of  u(x,t),  and  yields  the  amplitude  density  spectrum. 

A (W)  represents  the  amplitude  density  function  at  the 
source,  K the  attenuation  coeficient  and  0 the  epicentral 
distance  degrees. 

In  this  form  (Ewing  et  al.,  1957)  we  can  consider  a train 
of  travelling  surface  waves  as  formed  by  a synthesis  of 
single  harmonic  waves.  Let  the  suscripi;  1 refer  to  some  re 
ference  distance,  then 

\l(xifU))-  AM  . *xp  [-  K<4  j 
Vsin  ©, 

and 

■ 

V(*,w)  _ eXp[-K(x-x,)] 

\1  (x„fci)  \Z«h®/sin©-, 
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Talcing  the  logarithm  of  both  sides  of  the  preceding  equa- 
tion, gives 


V(*,u>)  -t n\/(xn4/)  = -i  /nSm6(  - A-fnSinfr-K'D 


where 


D=  x-x, 


Equation  (2)  can  be  written  as 


2 fix  V(x,4>)  + fn  sin  $ - C-  2 KD 


For  every  given  value  of  W,  the  terms  on  the  left  of  equa- 
tion (3)  contain  the  amplitude  density  at  a distance  X,  co 
rrected  for  the  effect  of  geometrical  spreading  over  the 
surface  of  the  spherical  Earth.  If  the  sum  of  the  termg  on 
the  left  is  plotted  against  the  independent  variable  D,  the 
resulting  straight  line  will  have  a slope  equal  to  -2K  and 
intercept  of  C.  The  equation  (3)  thus  provides  a method  of 
determining  K as  a function  of  frequency 


DATA  ACQUISITION 


By  visual  inspection  of  the  high  gain  photographic  records 
from  the  stations  used  in  this  paper  we  searched  digital 
tapes  from  1971  to  1974  and  selected  those  with  a clear 
train  of  surface  waves.  In  this  selection  we  eliminated  the 
records  with  significant  noise  in  the  frequency  range  of 
our  analysis,  i.e.,  0.01  - 0.1  Hz. 

Once  the  data  was  preselected  a computer  program  was  used 
to  get  the  digital  data  out  from  the  tape.  Again,  a selec- 
tion of  those  records  without  spurious  spikes  gave  us  the 
^mal  data  to  be  us«»d  in  the  present  study.  In  this  stage 
we  used  the  programs  HIGAN  and  EDITR  from  Lamont-Doherty  Ob 
servatory  to  obtain  the  data  off  the  tape. 


In  the  case  of  WWSSN  records,  the  normal  procedure  of  dig^ 
tizing  70  mm  films  provided  by  the  NOAA,  was  used  with  a di 
gitizing  interval  of  1 second. 

To  verify  the  digital  data  both  in  the  Standard  and  HGLP  re 
cords  we  developed  a computer  program  to  present  to  data 
with  the  computer  printer. 

We  were  able  to  check  the  errors  introduced  by  manual  digi- 
tizacidn  on  the  Standard  records  or  in  automatically  digiti 
zed  data  of  the  tapes  for  the  HGLP  case. 


COMPUTATION  OF  THE  AMPLITUDE  SPECTRA 

A computer  program  for  the  calculation  of  the  Fourier  trans 
form  was  adapted  to  obtain  the  data  from  the  tapes  in  the 
case  of  the  HGLP.  Again,  the  program  of  plotting  the  spec- 
trum was  used  but  for  a better  presentation  we  redraw  the 
points  joined  by  straight  lines.  Examples  of  such  spectra 
are  shown  in  figures  1 and  2 for  the  Standard  records  and 
figures  3»  4-  and  5 for  the  HGLP  system. 

After  the  Fourier  amplitude  spectrum  of  the  seismogram  has 
been  computedt  must  be  corrected  for  the  instrumental  ef- 
fects if  the  stations  have  instruments  of  different  charac- 
teristics. To  determine  the  attenuation  coefficient  K it  is 
not  necessary  to  correct  for  instrumental  effects  if  the 
seismographs  at  all  stations  have  the  same  characteristics, 
or  if  the  method  of  observaing  several  succesive  passages 
of  the  same  wave  at  a single  station  is  used.  We  will  dis- 
cuss later  on  about  this  method  which  clearly  point  to  the 
HGLP  system  as  advantageous  if  an  appropiate  velocity  fil- 
ter is  used.  In  the  present  study,  instrumental  corrections 
are  necessary,  because  the  attenuation  is  calculated  bet- 
ween some  reference  distance  and  the  observing  station  and 
also  because  all  seismographs  do  not  have  the  same  respon- 


ae  characteristic,,.  The  latter  point  is  especially  impor- 
tant when  seismograms  of  both  the  TOSS  and  Canadian  net- 
works are  used  or  HGLP  from  different  stations,  ^the  case 
of  WWSS  stations,  due  to  the  homogeneity  of  the  system  and 
similarity  in  the  shape  of  the  system  response  curves,  the 
magnifications  differ  only  by  a multiplier  which  does  not 
depend  upon  the  period  except  for  the  systems  which  have  a 
peak  magnification  at  6000. 

For  the  Canadian  network  and  HGLP  systems  the  corrections 
are  more  complicated  due  to  the  different  response  of  each 
instrument.  The  instrumental  response  of  each  station  was  ob 
tained  from  the  Bulletin  published  by  the  Dominion  Observa-"" 
tory,  Ottawa  and  the  list  of  amplifications  cf  the  HGLP  of 
the  Lamont-Doherty  Obe'’~vatory,  New  York. 


ATTENUATION  COEFFICIENT  AND  AMPLITUDE  SPECTRA  AT  a REFEREN- 
CE DISTANCE  “ ~ 


Following  the  approach  previously  described,  equation  (3) 
is  used  to  calculates  K.  A program  for  least-squares  appro- 
ximation of  a function  by  a polynomial  of  up  to  12  degrees 
was  adopted  for  a linear  fitting,  linked  with  the  main  pro- 
gram to  obtain  values  of  the  attenuation  coefficient  K,  the 
constant  C and  the  standard  deviation  of  K.  This  fit  is  ob- 
tained for  each  frequency. 

Due  to  some  window  effects  in  the  Fourier  analysis  program 
and  the  frequency  response  of  the  seismograph,  the  amplitu- 
des to  be  fitted  at  both  edges  of  the  selected  frequency 
band  are  not  very  reliable.  Consecuently  poor  accuracy  may 
be  obtained  for  K and  C at  those  frequencies. 

The  values  of  K,  plus  or  minus  one  standard  deviation,  are 
plotted  in  figures  6 and  7 to  give'  an  idea  of  the  range  of 
probable  values  of  K. 


In  equation  (3),  the  constant  C is  the  intercept  of  the 
straight-line  fit  to  the  amplitude  spectrum  at  a frequen 
cy  W at  some  reference  distance  when  a correction  is~ 
applied  for  geometric  spreading  from  the  source  to  x , 
uut  no  reduction  is  made  for  anelastic  attenuation  from 

the  source  to  A particular  value  for  x does  not  need 
to  be  specified. 

A plot  of  such  spectrum  provides  a possible  method  to  dis 
criminate  nuclear  explosions  from  earthquakes  as  Tsai  and 
Aki  (1971  have  shown.  There  is  a marked  difference  bet- 
ween the  spectra  corresponding  to  both  sources. 


ANALYSIS  OP  RESULTS  AND  CONCLUSIONS 

The  proposed  method  for  calculating  attenuation  coeffi- 
cient requires  a good  azimuthal  coverage  of  stations. 
This  condition  is  satisfyed  by  the  WWSSN  but  not  so  well 
by  the  HGLP;  the  problem  is  further  complicated  by  the 
failure  of  operation  of  some  stations  at  the  time  of  the 
explosion  or  by  the  fact  that  the  calibration  curve  was 
not  properly  known  or  very  poorly  defined. 

In  Figure  8 shows  the  average  value  of  the  attenuation 
coefficient  K for  Europe,  Canada  and  the  Eastern  United 
states  as  a function  of  frequency  from  the  1966  ,and  1970 
explosions  and  Standard  data.  A mean  value  of  K of 
loO  x 10  Km  is  found  in  the  frequency  range  of  O.Ol 
Hz  to  0.045  Hz.  Prom  here  there  is  an  almost  linear  in- 
crease up  to  0,10  Hz  and  about  0.075  - 0.080  Hz.  A mini- 
mun  at  0.05  to  0.06  Hz  .has  been  observed  previously  by 
Tryggvason  (1965)  and  Tsai  and  Aki  (1969)  and  a minimum 
at  about  0.08  Hz  is  consistent  with  that  also  found  by 
Tryggvason  (1965).  This  latter  minimum  in  the  data,  as 
observed  in  the  present  paper,  results  specifically  from 
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a good  transmission  of  energy  to  Canada  and  the  Eastern  Uni- 
ted States  in  the  frequency  range. 

It  is  observed  that  the  polynomial  fit  is  very  sensitive  to 
the  range  of  distances  between  stations,  Mezcua  (1973).  In 
this  sense  we  observe  less  scatter  in  the  results  whenVr- 
icing  with  stations  spread  out  over  large  distances  rather  than 
the  same  number  of  stations  being  concentrated  in  a small  dis 
tance  interval.  In  the  case  of  HGLP  data  a few  atations  were 
used  in  the  fitting  for  K calculation,  giving  values  with  a 
large  standard  deviation.  For  this  reason  we  decide  not  inclu- 
de them  with  the  set  of  the  data  because  the  lack  of  homoge- 
neity. An  alternative  to  use  the  HGLP  data  for  this  pourpose 
is  to  determine  K values  by  comparing  succesive  pass  throught 
the  same  station.  With  this  method,  used  by  many  authors  with 
Standard  Stations  it  was  necessary  to  use  large  magnitude 
earthquakes  and  no  nuclear  ezplosions.  Among  other  studies 
with  this  method  are  Anderson  at  al  (1965),  Bath  and  L6pez 
Arroyo  (1962),  Ben-Menahem  (1965),  Brune  (1962),  Oliver  (1962), 
etc.  Now,  with  the  HGLP  data,  however  will  be  possible  to  get 
down  with  the  magnitudes  to  be  in  the  explosion  range,  enough 
to  investigate. 

In  relation  with  the  amplitude  spectrum  at  a reference  distan 
ye»  the  use  of  Standard  Stations  give  us  a good  agreement  bet 
ween  the  theoretical  curve  (solid  line)  and  the  experimental 
(dotted  line)  in  Figure  9.  This  curve  corresponds  to  the  1970 
explosion*  In  this  sense  the  HGLP  data  shows  a clear  similar! 
ty  with  the  one  gathered  with  the  Standard  Stations. 
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TABLE  1 


LIST  OP  EXPLOSIONS  USED  WITH  HGLP  STATIONS 


DATE 

ORIGIN  TIME 

COORDINATES 

MAGNITUDE 

27  September  1971 

05 

59  55.2 

7i.4°N 

55.1°E 

6.4 

28  August  1972 

05 

59  56.5 

73.3°N 

55.1°E 

6.3 

21  September  1972 

15 

30  00.2 

37.08°N 

1 1 6 . 037°W  5.7 

12  September  1973 

06 

59  54.3 

73.3°N 

55.2°E 

6.8 

27  October  1973 

06 

59  57.4 

70.8°N 

54.2°E 

6.9 

LIST  OP  EXPLOSIONS  USED  WITH  STANDARD  STATIONS 


DATE 


ORIGIN  TIME  COORDINATES  MAGNITUDE  (MS) 

I 

27  October  1966  05  57  57.7  73.4°N  54.9°E  6.3 

H October  1970  05  59  57.1  73.3°N  55.1°E  6.7 
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- LIST  OP  FIGURES  - 


PIG.  1. 
PIG.  2. 
PIG.  3. 
PIG.  4. 
PIG.  $. 

PIG.  6. 
PIG.  7. 

PIG.  8. 
PIG.  9. 


- Amplitude  spectra  for  the  vertical  component  at  FLO 
(WWSS)  of  a Nuclear  Explosion  at  Novaya  Zemlya  (1970) 

- Amplitude  spectra  for  the  vertical  component  at  ALE 
(WWSS)  of  a Nuclear  Explosion  at  Novaya  Zemlya  (1966), 

- Amplitude  spectra  for  the  vertical  component  at  CHG 
(HGLP)  of  a Nuclear  Explosion  at  Novaya  Zemlya  ( 1 973 ) . 

- Amplitude  spectra  for  the  vertical  component  at  KIP 
(HGLP)  of  a Nuclear  Explosion  at  Novaya  Zemlya  (1972), 

Amplitude  spectra  for  the  vertical  component  at  OGD 
(HGLP)  of  a Nuclear  Explosion  at  Novaya  Zemlya  ( 1 971 ) 
incorrected  for  the  instrument  response. 

Values  of  attenuation  coefficient  for  the  1966  and 
197<J  explosions. 

Values  of  attenuation  coefficient  as  determined  by 
stations  in  Canada  and  the  Eastern  United  States  for 
the  1970  explosion. 

- Amplitude  spectrum  at  some  reference  distance  for  Oc- 
tober 27,  1966. 

- Amplitudespectrum  at  some  reference  distance  for  Octo 
ber  14,  1970 
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DEEP  STRUCTURE  OF  THE  MEDITERRANEAN  SEA  BY  SURFACE  WAVES  DISPERSION 


Gonzalo  Payo,  Elisao  Ruiz  da  la  Parte 
end  Jesus  Sierra.  Obaervatorio  de  Toledo. 


INTRODUCTION 


Though  in  the  last  years  several  papers  have  been  written  about  the 
nature  of  the  crust  and  upper  mantle  of  the  Mediterranean  Sea,  the  structure 
of  this  area  is  still  not  well  known  in  detail.  Figure  1 summarizes  most  of 
the  works  concerning  the  upper  structure  of  the  Mediterranean  Sea,  showing 
the  models  proposed,  by  different  authors  (see  References),  for  the  Medite- 
rranean crust,  according  their  results  obtained  by  surface  wave  studies  or 
explosion  seismology.  Lately  some  papers  have  been  also  made,  summarizing 
the  knowledge  about  the  whole  structure  of  the  region  and  its  geologic  evo- 
lution (Pyon  et  al.  1969j  Finetti&  Morelli  1973,  Payo  1975),  The  mediterra- 
nean area  appears  very  complicated  and  it  cannot  be  considered  as  a homoge- 
neous structure  since  the  inner  seas  (Balearic,  Alboran,  Tirrene,  Jonic,  Adria 
tic,  etc)  have  had  different  origins  and,  presently,  their  crustal  structure 
shows  also  notable  differences.  In  some  of  our  previous  papers  (Payo  1967, 
1969)  we  found  a great  structural  difference  between  the  West  Mediterranean, 
with  a crust  almost  oceanic,  and  the  East  Mediterranean  which  appears  with 
a crust  of  a more  continental  type,  having  a thick  sedimentary  layer. 

One  of  the  parts  less  known  of  the  Mediterranean  structure  and 
possibly  the  most  homogeneous  one,  is  that  corresponding  to  the  upper  mantle. 
That  is  the  reason  why  the  present  work  try  to  study  this  deep  structure 
along  some  mediterranean  paths,  by  using  the  surface  wave  dispersion  methods. 
The  Rayleigh  wave  phase  velocity  will  he  studied,  considering  the  propagation 
of  the  surface  waves  between  each  pair  of  long  period  stations  located  in  the 
Mediterranean  borders. 


In  the  first  part  of  the  work,  we  have  used  the  two  high  gain  long 
period  stations  (HGLP)  existing  in  the  region,  Toledo  (Spain)  and  Eilat 
(Israel).  The  characteristics  of  these  stations,  of  a very  long  period  and 
high  gain  (Savino,  J.M.  et  al.  1972)  offer  important  advantages  for  this 
kind  of  studies.  First  of  all  their  particular  location  permits  to  study  a 
long  profile  (TLO-EIL)  that  crosses  longitudinally  the  whole  Mediterranean 
area  (see  Figure  l).  The  long  period  of  this  system  (Tg=100  sec.,  Tp^O  ~cc. ) 
makes  possible  to  record  Rayleigh  waves  of  a long  period,  even  from  earth- 
quakes of  moderate  magnitude.  Finally,  the  agreement  of  the  path  TLO-EIL  with 
the  ones  between  the  WWSS  standard  stations  TOL-HLW  and  TOL-JER , permits  to 
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compare  the  information  given  by  both  system  of  equipments  (HGLP  and  WWSS)  and 
discuss  the  advantages  relative  to  each  of  them,  for  the  study  of  the  crust- 
mantle  structure. 

A second  part  of  this  work  was  scheduled  to  study,  the  other  mediterra- 
nean paths  that  can  be  formed  by  using  pairs  of  standard  stations  located  on 

the  borders  of  the  Sea  (as  Toledo,  Malaga,  Aquila,  Jerusalem,  Helwan  and 
Athens] . 

In  this  report,  however,  only  the  part  of  the  results  concerning  to  the 
HGLP  stations  is  shown,  because  the  data  of  these  station  are  the  only  ones  we 
have  had  time  enough  to  process.  Nevertheless,  the  reproduction  and  digitation 
□f  the  standard  seismograms  has  been  already  made  for  all  the  pairs  of  standard 

stations.  This  report,  therefore,  only  shows  some  provisional  result  of  the 
whole  work. 

We  must  say  that  the  rhythm  of  the  work  has  been  interrupted  several  times 
by  different  causes.  For  instance,  the  difficulty  of  getting  acceptable  data 
from  the  digitation  of  the  magnetic  tapes,  filed  in  the  Lamont-Doherty  Obser- 
vatory, where  we  went  to  pick  up  them.  This  work  was  scheduled  from  the  begin- 
ning, by  assuming  that  the  magnetic  tapes  of  the  HGLP  system  cculd  be  readily 
used  to  get  the  digitation  of  the  wave  trains.  The  impossibity  of  using  the 
tape  data  has  led  to  the  necessity  of  digitizing  afterwards  the  seismograms 
directly  from  enlarged  microfilm  copies,  with  the  corresponding  delay  in 
the  work  to  doing  that.  Other  additional  trouble  for  this  cause,  has  arised 
with  the  difficulty  in  digitizing  large  amplitude  seismograms,  which  had  been 
avoided  if  the  digitation  had  been  obtained  from  the  magnetic  tapes.  The  exis- 
tence of  these  large  amplitudes  has  forced,  sometimes,  to  reject  the  fundamen- 
tal mode  in  several  cases.  Another  delay  in  the  work  was  produced  by  the 
change  of  computer  in  the  Center  were  the  programs  were  being  processed,  since 
the  IBM  7090  was  substituted  by  an  IBM  360/65.  Therefore,  when  the  programs 
had  been  prepared  for  the  first  computer  we  had  to  repeat  all  the  process  of 
adaptation  of  the  programs  to  the  new  computer.  That  substitution  produced 
also  a period  of  more  than  three  months  of  inactivity  in  the  Processing  Center. 


THE  DATA 

All  the  earthquakes,  whose  epicenters  were  lined  up  with  the  stations 
TLO  and  EIL  have  been  examined  from  1971  (date  of  installation  of  the  HGLP 
system]  up  to  the  present.  Only  six  earthquakes  for  the  profile  Toledo-Eilat 
and  three  for  the  reverse  profile  were  finally  accepted.  The  optimum  records 
for  this  kind  of  studies,  where  the  Rayleigh  fundamental  mode  is  used,  corres^ 
pond  to  earthquakes  which  epicentral  distances  are  about  5 to  10  thousand  km 
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from  the  closest  station.  In  this  sense  the  seismic  zones  linen  up  with  Toledo 
and  Eilat  (South  America  and  the  Indian  Ocean)  are  not  seismically  very  active 

(specially  the  last  one),  which  makes  difficult  to  find  a large  number  of  appro- 
pi  iate  seismograms. 


In  most  cases  the  Rayleigh  fundamental  mode  has  been  used,  being  the  two 
stations  considered  at  some  15  wave  lengths  of  distance,  relative  to  the  largest 
period  visible  on  the  seismogram.  This  distance  is  somewhat  large  with  relation 
to  the  wave-length,  since  the  optimun  values  correspond  to  distances  about  3 to 
6 wave-lengths.  However,  the  filtering  techniques  and  the  methods  of  cross- 
multiplication for  determining  the  phase  velocity,  may  Dead  to  good  results. 

Also  we  have  used  waves  of  great  periods,  as  the  surface  waves  R2  and  R3  (Ray- 
leigh wave  giving  turns  around  the  Earth),  to  compensate  the  large  distance 
between  Toledo  and  Eilat. 


Most  of  the  digitation  of  the  seismograms  has  been  made  with  an  automatic 
digitizer,  that  takes  the  readings  from  the  enlarged  microfilm  copies,  when 
we  follow  the  signal  by  means  of  an  optic  pointer.  The  readings  have  been  made 
every  second  for  the  fundamental  mode  and  every  two  seconds  for  R2  and  R3  waves. 
The  pairs  of  seismograms  corresponding  to  the  standard  stations  (13  cases  in 
total)  have  been  digitized  in  the  same  manner. 


From  the  nine  earthquakes,  above  mentioned,  corresponding  to  the  HGLP 

stations  (TLO-EIL),  four  or  them  shown  R2  and  R3  wave,  besides  the  fundamental 

mode.  However,  some  of  these  cases  were  unfortunately  rejected,  due  to  the 

difficulty  in  digitizing  the  large  amplitudes  of  some  of  the  seismograms. 

Therefore,  the  subsequent  analysis  has  been  only  applied  to  five  cases  of 

fundamental  mode,  five  of  R2  waves  and  one  case  of  R waves. 

o 

A program  for  multiple  filtering  (Dziewonski  et  al.  1969)  was  applied 
to  all  these  cases.  The  program  reduces  the  digitized  values  to  the  average 
of  them,  removes  the  linear  trend,  computes  the  Fourier  transform  (algorithem 
of  Cooly&Tukey)  and  applies,  along  the  seismogram,  a gaussian  filter  truncated 
at  30  db  and  with  a bandwidth  proportional  to  the  frequency.  Finally,  the 
program  gives  the  energy-density  diagrams  as  a function  of  the  group  velocity 
and  the  logarithm  of  the  period.  Figures  3a  a 3b  show  an  example  of  the 
diagrams  obtained,  corresponding  to  the  same  earthquake  recorded  at  the  two 
stations. 


Regarding  the  results  of  all  these  group  velocity  diagrams,  some  of  the 
cases  with  little  energy  or  with  the  contour  lines  not  well  defined  have  been 
rejected.  The  rest  of  the  cases  accepted  appear  in  the  next  table. 


TABLE  I 


DATE 

ORIGIN 

TIME 

h m s 

x 

1 

FIRST 

STATION 

SEGONO 

STATION 

TYPE 

OF 

WAVE 

\ 

REGION 

7-Sep-72 

□2-54-58 

22O&-602QE 

EIL 

i 

TLO 

i 

F.  mode 

Indian  Ocean 

16-Sep-72 

09-14-34 

159 2N— 9622W 

TLO 

EIL 

F.  mode 

Mexico 

23-Dec-72 

06-29-42 

1294N-B691W 

TLO 

EIL 

F.  mode 

Managua 

7-Jun-73 

1&— 34— 46 

1492N— 9199W 

TLO 

EIL 

F . mo  de 

Guatemala 

lG-Sep-72 

09-14-34 

1S92N-9692W 

EIL 

TOL 

r2 

Mexico 

14-Apr-73 

08-34-00 

1097N-8427W 

EIL  | 

TOL 

r2 

Costa  Rica 

Ths  group  velocity  curves  thus  obtained  for  all  these  cases  have  been 
used  to  get  the  values  of  the  phase  velocity  across  the  path  between  Eilat 
and  Toledo.  For  the  computation  we  have  used  a program  of  cross-multiplication, 
FACRQS,  that  from  the  filtered  seismograms  at  the  two  stations  and  by  using 
the  group  velocity  previously  obtained,  computes  the  phase  velocity  for  the 
range  of  periods  14  to  140  seconds.  In  general  lines,  the  program  operates  as 
follows.  First,  it  computes  the  arrival  time  of  each  period  for  both  seismograms, 
in  order  to  apply  a "window"  with  a width  five  times  the  corresponding  period. 
Next,  for  each  one  of  the  vel  cities  between  3 and  5 km/sec,  the  program  com- 
putes the  time  = (A-B)/cj_  and  it  moves  one  seismogram  relative  to  the 
other  the  amount  T±,  making  in  this  position  the  multiplication  of  them.  The 
largest  of  these  products  is  plotted  against  the  period  and  phase  velocity  in 
form  of  a normalized  diagram  (see  the  example  of  figure  4).  The  phase  velocity 
curve  appears  then  as  the  line  that  joints  the  maximum  values  of  the  diagram. 

We  have  designed  with  thick  circles  the  highest  values  of  each  column  and  with 
thin  circles  the  other  large  values  close  to  the  maximum  . 

RESULTS  AND  DISCUSSION 

The  phase  velocity  curves  for  all  the  cases  finally  admitted  have  been 
averaged.  Figure  5 shows  this  average  (doted  line),  where  the  vertical  bars 
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represent  the  standard  deviation  for  some  selected  periods.  In  the  figure  we 
have  also  plotted  the  curves  corresponding  to  several  models  to  be  compared 
with  our  results. 

As  we  can  see,  the  curve  here  obtained  shows  fairly  low  velocitfesfor 
the  range  of  periods  between  15  and  35  seconds.  This  part  of  the  curve  is  the 
most  reliable  one,  because  the  velocities  are  in  good  agreement  for  all  the 
cases.  These  phase  velocities  are  governed  by  the  shear  velocities  of  the 
first  layers  of  the  Earth,  up  to  some  60  kms  of  depth,  and  these  low  values 
were  already  found  by  the  author  (Model  EME  of  the  Figure  ) for  the  eastern 
Mediterranean  Sea.  Of  course,  the  present  values  are  somewhat  higher  than 
those  of  model  EME,  but  we  have  to  consider  that  the  path  TLO-EIL  only  has  a 
SO  per  cent  along  this  eastern  part  of  the  Mediterranean.  However  the  influence 
of  this  part,  with  this  thick  sedimentary  layer  that  seems  to  have  (Payo  1969), 
is  clearly  visible  in  the  dispersion  curve,  whose  values  in  the  short  period 
range  are  below  those  corresponding  to  a typical  oceanic  model  (8099  in  the 
Figure)  or  a continental  model  (CANSD). 

The  rest  of  the  dispersion  curve,  here  obtained,  that  is  the  part  of 
the  curve  between  40  and  110  seconds  of  period  exhibits,  however,  fairly 
high  values.  The  phase  velocities  are  comparable  to  those  of  the  CANSD  model, 
and  they  are  some  0.1  km/sec  larger  than  those  of  the  oceanic  mudel  8099  and 
about  0.2  km/sec  larger  than  those  of  the  models  G1  or  G2  given  by  Berry  &. 
Knopoff  for  the  West  Mediterranean  basin.  The  interpretation  of  such  high 
velocities  is  problematic.  It  would  suppose  high  shear  velocities  below  the 
Mohorovicic  discontinuity  ( about  4.7  km/sec)  and  the  absence  or  little  signi- 
ficance of  the  low  velocity  channel  in  the  upper  mantle.  Other  posibility  is 
that  the  channel  be  located  at  a depth  greater  than  normal.  The  flat  nature 
of  the  curve  for  periods  larger  than  80  seconds  is  in  favor  of  this  hypothe- 
sis. 

These  first  results  will  be  complemented  with  future  studies  about  this 
zone,  and  they  will  give  new  information  about  the  structure  of  this  complex 
mediterranean  region. 

Gonzalo  Payo  Subiza 
Eliseo  Ruiz  de  la  Parte 

Jesus  Sierra  Abril,  1975. 

Observatorio  Geofisico  Central 
Instituto  Geografico  y Catastral 
TOLEDO 


FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4. 


Figure  5. 


Some  of  the  models  obtained  for  different  areas  of  the  Mediterra- 
nean sea.  LR  means  Rayleigh  waves,  LQ  Love  waves,  R refraction 
profiles,  G group  velocity  methods  and  F phase  velocity  methods. 

The  figures  without  parenthesis  are  the  velocities  given  directly 
by  the  method  of  investigation  used. 

Also  the  path  TOL-EIL,  studied  in  this  work,  appears  in  this  Fi- 
gure. 

i-  An  example  of  the  seismograms  used.  It  corresponds  to  the  earth- 
quake of  7 Sept.  1972,  recorded  at  TLO  and  EIL  with  the  HGLP  seis- 
mometers. 

,-  Group  velocity  curve  obtained  "multiple  filtering"  method.  The 
curve  appears  as  the  line  that  follows  the  maximum  energy  of  the 
diagram.  It  corresponds  to  the  earthquake  of  16  Sept.  1972,  recorded 
at  the  TOL  (Figure  3a)  and  EIL  (Figure  3b). 

- Diagram  of  the  phase  velocity  curve  between  TOL  and  EIL  obtained  by 
the  cross-multiplication  method  (see  the  text).  It  corresponds  to 
the  earthquake  of  7 Sept,  1972  in  Indian  Ocean  region. 

Comparison  of  the  phase  velocity  curves  here  obtained  with  those 
corresponding  to  some  Mediterranean  models  and  with  a typical  oceanic 
or  continental  structure.  Models  1 and  2 correspond  (Berry &Knopoff, 
1967)  to  the  West  Mediterranean  basin,  model  3 (Payo,  1969)  to  the 
Eastern  part  of  the  Sea,  model  4 is  an  Atlantic  model  (Dorman  et 
al.  I960)  and  model  5 a typical  continental  model,  (Brune  and 
Dorman  1963), 
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